Non-Gaussian states produced by close-to-threshold optical parametric oscillators: 

role of classical and quantum fluctuations 
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Quantum states with non-Gaussian statistics generated by optical parametric oscillators (0P0) 
with fluctuating parameters are studied by means of the Kurtosis excess of the external field quadra- 
tures. The field generated is viewed as the response of a nonlinear device to the fluctuations of laser 
pump amplitude and frequency, crystal temperature and cavity detuning, in addition to quantum 
noise sources. The Kurtosis excess has been evaluated perturbatively up to the third order in the 
strength of the crystal nonlinear coupling factor and the second order in the classical fluctuating 
parameters. Applied to the device described in Opt. Expr. 13, 948-956 (2005) the model has given 
values of the Kurtosi excess in good agreement with the measured ones. 
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I. INTRODUCTION 



Non Gaussian (NG) resources (i.e. NG states and/or 
operations) [3-@] are required for realizing relevant quan- 
tum information protocols, as for example entanglement 
distillation and swapping (7rlToj. It has been demon- 
strated that they may improve the fidelity of teleporta- 
tion [TlH3 and cloning [3] and NG states are more ef- 
fective in revealing nonlocality (l5rll8j . Thus the reliable 
generation of NG single- and two-mode states assumes 
a relevant role. Recently, few scheme to generate NG 
states based on conditional de-Gaussification protocols 
have been proposed [liT - [ill ] and realized [l^, H(| ■ More- 
over, it has been proven that phase diffusing a squeezed 
vacuum state makes it a NG ones 



Departures from Gaussian statistics have been ob- 
served (see Ref. [2l|) in the outcomes for the field quadra- 



ture Xg = (e 



/2 outing a degenerate Opti- 



cal Parametric Oscillator (OPO). These deviations have 
been quantified by measuring for different operating con- 
ditions the departure of the fourth moment (A 4 ) from 
its Gaussian value. In particular, it has been observed 
a dependence on the quadrature phase 9 with the max- 
imum departure always appearing in correspondence of 
the anti-squeezed amplitude quadrature (0 = 0). 

No doubt the observed deviations are due to a non- 
linear response of the device to Gaussian quantum and 
classical fluctuating parameters, such as laser pump am- 
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plitude and phase, cavity length and OPO crystal tem- 
perature. While in the linear analysis only the quantum 
input noise contributes to the OPO output signal, by ex- 
panding the equation of motion to higher orders classical 
fluctuations contribute as well to the output. As a result 
of the multiplicative mixing of several noise sources, the 
signal looses its Gaussian character. The weight of these 
sources in determining the deviation from the Gaussian 
statistics was not clear at time of publication of Ref. [2l[ . 
The present paper is meant to present and discuss a con- 
sistent theoretical framework to explain and interpret the 
observed experimental evidences. 

OPOs rely on parametric down-conversion: a strong 
pump beam at frequency u> p interacts in a non-linear 
crystal with the vacuum fields thus generating two 
beams, signal and idler, at frequencies u> s and u>i re- 
spectively [22l - |24j ]. This mechanism is represented in 
the Hamiltonian by the product of three field operators, 
a p (pump), a s (signal) di (idler) and described dynam- 
ically by Langevin equations. The model herein pre- 
sented starts by including in the Graham and Haken 
Langevin Equations (GHLE) different classical noise 
sources. Then, expanding a s ^ as power series in the 
strengths of the quantum and classical fluctuations a 
hierarchy of Langevin equations is obtained in which 
the field at a given order acts as source for the next 
one. This has been done up to the third order thus 
culminating with Fig. [5] showing an adequate agree- 
ment of the computed difference (Kurtosis excess) Kg = 

((Xf) - 3 (Af ) 2 ) / (X|) 2 with the measured values K e 

of 

A perturbative analysis of an ideal OPO was already 
developed in (26j : according to it nonlinear contributions 
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become comparable to the linear output only at a rela- 
tive distance from the threshold of about ~ 10 ~ 6 while 
the data of Ref. [2l[ were measured at a distance of 
~5x 10 -2 so signalling the occurrence of more complex 
mechanisms. Hence, it has been essential to account for 
contributions from different fluctuating parameters each 
represented as a Gaussian process weighted by its stan- 
dard deviations gi. For the sake of generality the GHLE 
[iU model has been developed for a non degenerate OPO 
whereas the field moments have been calculated for a de- 
generate one so to allow a direct comparison with [2l[ . 

Different order fields have been represented as the con- 
volution of the previous order one with a 2 x 2 matrix 
G which becomes singular on approaching the threshold. 
This singularity is at the origin of the enhancement of 
higher order effects in the critical region. 

The paper is organized as follows. In Section [TT] the 
extended GHLE model for an OPO with fluctuating pa- 
rameters is introduced and discussed. The squeezing of 
the intracavity field is discussed in Sec. IIII1 Section [TV] 
deals with the statistics of the quadratures outside the 
cavity as measured by a finite bandwidth detector. The 
deviation from the Gaussian statistics in terms of Kg 
is analyzed by dwelling on the agreement of the results 
provided by the model with the experimental findings. 
Plots of Kg vs. 6 and Kg for 9 = vs. E 2 are reported. 
Section [V] closes the paper with concluding remarks. De- 
tails about the linearization, together with few analytical 
derivations are reported in the appendices. 



II. LANGEVIN EQUATIONS 

Consider the set of the three OPO cavity modes a k 
(a p = e~^ p * _l * p ao pump mode at frequency io p and 



(GHLE) gjj: 



phase 



iuj s t- 



*5'P» , o 1 and di = e 



respectively signal and idler modes with uj p = u) s + Wj) 
whose mutual interaction under the action of a driving 
field e~ lUpt £ is described by the Hamiltonian 



Hi, 



ih 



X 



0*10,2% 



aoa\a2 



ih (, 



£*a 



(1) 



where x is the coupling parameter proportional to the 
crystal second order susceptibility x^ ■ Since real lasers 
are characterized by a field £ = e {l + g Pp fi p ) e~ l< ^ p 
of constant amplitude e modulated by a fluctuating 



1 

-i<t>p 



SmpMpW ((fh) = 
where <p p (t) is a 



0) times a phase fac- 
slowly diffusing phase, 



4>v (*'))) = A e \t-t'\ with A e the laser 



factor 
tor e" 

i.e. (0 P (£) 
linewidth|24 

The cavity modes are characterized by damping factors 
7fc,M, Jk,x ilk = 7fc,M + 7fc,x) due respectively to the out- 
put mirror (M) and the other loss mechanisms (x: crystal 
absorption and scattering, absorption of the two mirrors, 
etc.). The evolution of the cavity mode operators can 
be described by the Graham-Haken Langevin equations 



Id ~ w i - 1 



, * t 



j = 1,2 and j ^ f 



a 



vy - ivo - i(j) p ) a - -y a i a 2 + e (l + .g Mp M P ) 



(2) 



where R k (t) = ^J2^ kM b kM + ^J2^ k}X b Kx takes 
into account the delta correlated vacuum fluctuations, 

b k ,M,x (*) (O) = Ht-t'), entering the OPO 

cavity. Modes are assumed to be slightly detuned by 
v k = f£ [L k ^\ - u k , with L k the OPO optical length 
at frequency ui k and [x] the closest integer to x. 

In the following we will indicate by n k = "f k — iv k = 
\n k \e~ v>f ' k (with ipi = ip2) the complex damping coeffi- 
cients, by k = \ki + K2I /2 the mean decay rate and by 
t = Kt the time normalized to the cavity lifetime k . A 
caret will mark quantities normalized to k (e.g. e = e/«) 
and a tilde those such that the integral of their correlation 
function (e.g. (£l p (t) fi p (r')) = C Mp (r - t')) is equal to 
1 (e.g. C Pp (r) dr = 1). In particular the Gaus- 

sian delta-correlated process d(f> p /dt will be replaced by 
d4> P /dr = g^ p w p with (zu p (t) zu p (r')) — S (t — r') and 



In real devices, beside the fluctuations related to clas- 
sical noise of the laser beam, the parameters v k and x 
of Eqs. ([2]) experience also the effects of mechanical 
vibrations. Residual fluctuations of the cavity optical 
length SL k , at frequency u k , induce deviations 5v k — 
— (5L k /L k ) uj k of the mode detunings from their average 
values (i> k ) . Usually, an active control guarantees that 
the standard deviation of Svq is a small fraction of 70- 

The parameter x 18 proportional to the crystal sus- 
ceptibility x^ through the Boyd-Kogelnik function 
Hbk (&, x, £,) of the phase-matching factor a, the fo- 
cusing parameter £ and the crystal absorption h. a (T) 
depends on the crystal temperature T through the refrac- 
tive indices at the interaction wavelengths. If the cavity 
configuration is far from the concentric one the depen- 
dence of £ on the cavity geometry fluctuations can be 
neglected. Under this assumption x will be replaced in 

the system ^ by xe _l( ^ x fl + 9t5T^J with x depending 
on the slow variations of T while <p x is a phase depend- 



ing on the position of the beam waist with respect to the 
crystal center. With an accurate alignment 4> x can be set 
equal to 0. gr is defined by: 



9t 



- d\ogH BK 
dT 



In conclusion, the OPO analyzed in the following is 
characterized by four classical fluctuating parameters 
5^ P Ap:5ra p ro P ,5^^ an d grST, where fi p , w p , 5v, and 
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ST are Gaussian processes with unit standard devia- 
tions, and g^g^, g Vk (= (g^/^) z%9» ) >9t the cor- 
responding weights. These four terms, together with 
g x = \\\/(2\k()k\) 1 ' 2 , describing the non-linear interac- 
tion of strength x [2a|: determine the OPO dynamics. 
For typical operating conditions, (g x ~ 10" 6 , K ~ 10-^20 
MHz, y/iAf) ~ 1 + 1000 Hz, x/iSf^) ~ 1 10 mK, 
and dn/dT ss 10" 6 10~ 4 ) g^ p , # rop , 5„ fc , g T range in the 
respective intervals 10~ 2 ~ 10" 1 , 10" 4 -f- 10~ 2 , 10~ 5 -j- 

io-\io- 5 -^io- 4 . 

The extended GHLE system @ may now be written 

as 

.9w p -& P . X A 
a j = - [Kj - i — - — + ig Vj 6v J a,j 

+ \[W( 1+9TSf ) 9 * a ° a i' 

+ e -ifl«„*p/2^. 

do = — (ko - igvj p ^ P + ig Vo Sv) ao 

- ( 4 + 9t8T) g x aia 2 + e (l + ff^Ap) 

+ e *fl-,*p£ (3) 

a dot indicating derivatives with respect to r. 

The OPO admits a threshold value for the amplitude 
e = e" 1 = |k | \/\kik 2 \/ (2 |x|). Below threshold, the 
mode ao has a non-zero mean value ro, which is related 
to the driving field amplitude e (e = «o^o)- Therefore 
separating the average part r from the fluctuating one 
Sao = ^oQ!o we put ao = ro (1 + ao), where ao = po — itpo- 
Conversely the modes a,j have zero mean value and will 
be expressed in terms of rescaled operators aj — rjaj, 
with the Tj defined in terms of ro 

\^\r] = \k \rl. (4) 

Passing now from the amplitudes at to the scaled quan- 
tities ctk, the extended GHLE © is rewritten as: 

dy + kjoij =E \kj\ a], + g x N x . 

+ \kj\E (gfST + a ) a], 

+ EgfSTaoOj, 
a + K a = - E |k | a x a 2 + 9 x ^xo + «o5^„Ap 

+ i9zu p ^p - ig Vo Su + igzup^pCto 

- ig„ Sva - E \h Q \ g T Sfa 1 a 2 , (5) 

where 

E= 9 \ = ee**» (6) 
V2k kik 2 

The above equations describe the dynamics of the fluctu- 
ating fields a^. as responses to the classical noise sources 



A P , vjp, 5u, ST and quantum terms 

N Xj =R j e- i ^*'/ 3 /(jg x r s ) 
N Xo =R o e^»/(g x r ) . 

The system © may be solved perturbatively upon ex- 
panding the field amplitudes as 

^ = E4 m) - (7) 

m— 1 

The first-order terms (m = 1) correspond to the lin- 
earized system. c*fc™~ 2 ' are generated by non linear 
sources sj." 1 ^ made of a quantum contribution, propor- 
tional to <7™, and of mixed terms involving products of 
quantum and classical fluctuations of the type g x _1 gi 

and g x ~ 2 gigj. Any involves fields calculated up to 
the (to — l)-th order. By substituting the a^ expansion 
into ([5]) and grouping the terms corresponding to the 

same perturbative order, each can be calculated as 
the convolution of the components Gkk' of the Green's 
matrix (see Eq. (|A1J> ) relative to the linearized system, 
with s[" 1 ' and sl™^. 

For a degenerate (the signal and idler fields collapse 
into a single field in this case) and tuned (vk — 0) OPO 
the extended GHLE (0 reduces to 

a + ha —Ea^ + g x N x + E (^gfST + a^j a^ 

+ * Qfcp^p - 9vj5P \ a + Eg f Sfa a^ 
a + h a =g x N Xo + £ g Mp Ap + ig™ P ™ P ~ igu Si> 

E 1 A 1 o • r» 

--^-|«o|a + «g rop ro p ao - i5„ <wa 
-£|£ 1 5TdTa 2 , (g) 

where the parameter 

E = -^=g = e*/ 2 4 = e^°/ 2 |£| 
V2^ I th 1 1 

now represents the excitation e normalized to the thresh- 
old e th = x/\k \ /2/g x while N x (t), N Xo (t) and w p (r) 
are delta-correlated processes, 

(iV x (r)iVt(r')) = ^(r-r') , 

(n xo (t)n1 (t>)\ = ^fs(r-T>) , (9) 

The correlation times for jx p are typically of the order of 
0.2-^1 /is, while those for Sv and ST are of the order of ms 
and s respectively. Although 8i> and 5T can be treated 
adiabatically, we have preferred to treat the noise sources 
in a unified fashion. 
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A. Nonlinear terms for a degenerate OPO 

Expanding Eqs. ([5]) up to the third order the fields 



(signal=idler) and a.^"^ (pump) of ([7]), repre- 



sented in the vector form a^— {o^ m \aS m ^Y ,ot^ = 



( (m) (m)t\ • i 

a (m /(*G£ - t) .^- m \ T )dT 



with G and Go defined by Eq. (| A2|) . while the signal 
s ( m ) anc j pump sources read respectively: 



S W(r)=.g x N x (r) 

s (2)( r ) =B W(r)-aW(r) 
s (3)( r ) =B (\(i)f( r ) 

2 



<?-<5£ (2) (r - t') 1 + B^ 1 ' 1 ) (r - r')J ■ a (1) (r') dr' 

(10) 



and 

4 1} ( T ) =5xN Xo (r)+i [Sw p ^ p (t) - 9uoSv(t)] 1_ 
+ «05mpAp(t) 

s o 2) ( T ) =i [ffn^ptto - ip,/ ^ao] 1- - --EkW 1)2 (11) 



with 1_ = (1,-1) T and BM 



B« - 



E* (Pa™ + g T 8f ) -i {\g^w p - gJD) 



(12) 

with P the permutation operator. Note that ST and ag 1 '' 
appear in the off-diagonal terms while w v and Sv are in 
the diagonal ones. 



J 



For a tuned OPO (V> = = X = °) fi(2) > <^ (2) and 
B^ 1 ' 1 ), of Eq. (JTUJ-c), are given by 



£( 2 ) = £ 



(2) 



2(1- i? 2 ) ya * £ 



2# 



^< 2 ) (r - /) = -S 2 « Go (r - r') (r - r') 
B (i,D ( T _ = / B (i) (r ) . G (r - r') • B« (r' 



with er£, 1) (T-T / ) = (aW (r) a™ ( r ')>. To be 
consistent with the above hierarchy, all moments 
(s( m V" l) ) must satisfy the inequalities (s (1) s (1) ) > 
(s (2) s (2) ) > (s( 3 )s( 3 )). Close to the critical point, 
(s^sW), (s( 2 M 2 )) and (s< 3 V 3 >) are of the order of 

O ( 5 2 ), O (glgj (1 - Ey 1 ) and O (firf (1 - £)- 2 ) re- 
spectively, so that the above condition implies that the 
approximation maintains its validity up to 1 — E > gf. 

Since the linear source for the down-converted beam 
(see Eq. (fTUl -a)) contains uniquely the quantum noise 
term, the effects of the classical fluctuations can be ana- 
lyzed only going beyond the linear approximation. 



III. INTRACAVITY FIELD 

The nonlinear contribution to the intracavity field is 
represented by the averaged tensor product with re- 
spect to the different fluctuating parameters g L N L {g x N x , 

9x N xo> 9» P fr P , 9ru P ^p, 9u5i>, and g T ST), 

(r) = (a (r) a T (0)) - (a« (r) a^ T (0)) 



0.(2,2) ( T ) +<T (3,1) ( T ) +0 -(1.3) ( r ) 



4. 



Relevant er^™ 1 -™) (r) are explicitly given by: 



r (2,2) 



(r) = 



f T dr 1 f dr"G (r - r') ■ (b« (r') • (r' - r") • B^ T (r")) ■ G (-r") 

J — oo J — oo 
r 



(r) = I drG (r - r') • (V> \ \ I } + f (<5B {2) (r - r') 1 + B^ 1 ' (r - r')) • <x (M) (r') rfr' 

J — oo \ L J J — oo 



Of particular interest is the variance {xfj^J nor- the nonlinear correction to the squeezed variance: 



malized to ( X 



L ir/2 



the weight, normalized to g 2 , of 



Af Vi = 



y(NL)2\ 



^/2 



<»'+^ t ) (Q)+< r'»'-<;:i'<») 

4^ 2 /2 



(N L) 



r (NL) 
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For a balanced and exactly tuned OPO we have 



A 



NL 

Xo 



A 



NL 



i NL 



v NL 



E (2E 3 k + 2E 2 k (2 + k ) + (2 + k ) 2 - 2E(-2 + k (2 + k ))) 
2(1 - E 2 )(l + £0(2 + k )(2 + 2E + k ) ' 
k Q {2 + k ) - Ek (6 + kp) + 2E 3 {8 + 5k Q ) - 2E 2 {\2 + k (9 + k )) 

16(1-£ 2 )/c (2 + k ) 
g (E 3 + 2k - £ 2 (3 + ftp) - £(6 + kq)) 
2(1-£ 2 )(1 + £)k 2 ' 
£ 2 (1 + EQ 2 - (2 + g) (1 - E)k 2 
2(1 - £ 2 ) (1 + E) k 2 
2 + E 2 
2{1+E) 2 ■ 



(13) 



As it is apparent from the above formulas the contribu- 
tion of ftp remains bounded on approaching the threshold 
[E — > 1), whereas the other terms diverge as (1 — E)^ 1 , 
see Fig. [TJ 
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FIG. 1: Nonlinear contributions A NL (see Eqs. {T3}) to intra- 
cavity squeezing, as functions of the normalized pump ampli- 
tude E for ko = 5 (left) and kq = 10 (right). In both plots: 
A^ L -> black solid, A" L -> black dashed, Aj? L -> black dotted, 



i NL 

A^ L — s> gray solid, and Ay L — > gray dotted. 



It it noteworthy that the ratio between A5£ and the 
analogous quantity A^ SE calculated in [26[ by means of 

A NL 

the Positive P representation (PPSE) goes as , P p° E ~ 

xo 

+ 0(1-E). In the limiting case of E ~ 1 the two 
approaches differ by a Ko-dependent factor bounded be- 
tween 1/3 and 1. Such a substantial agreement between 
the two results in proximity of the singular point vali- 
dates the use made in the present paper of the extended 
GHLE. 



IV. KURTOSIS EXCESS AND COMPARISON 
WITH THE EXPERIMENTAL RESULTS 

The held a out ,i outing the OPO is a function of a, 
the mirror damping coefficient 71 , and the corresponding 



input noise Ni |2£ 



OioutS 



27: 



(14) 



Accordingly the generic output quadrature reads 



Xg 




where 9 = | 



■ e -io e m 



x. 



) and ATg" 1 ' corresponds to a 



(1) 



X 



(2) 



^(™) 

While Xg 1 ^ is Gaussian the terms Xg m>1 ^ deviate from 
the normal distribution. 

Quadratures are detected by a balanced homodyne and 
the relative current is measured by selecting a frequency 
Slf and an integration time 1/7/ [2l],|30|. Accordingly, 
the detector output is represented by 



V e = F f X. 



2 ^ T ' cos (O/t') X g (r') dr' (15) 



The deviation of Vg from a Gaussian distribution can 
be measured by the Kurto sis- excess parameter Kg (see 
Appendix IDl for details) 



Kg = 



(:^:)-3(:^:) 2 
3(:V 2 :) 2 



(1)2 



(16) 



with Tg L given in Eq. (|Dip . the weight of the different 
noise sources indicated generally by N L . 

The spectral density S^ p (w) = (n^p (—to) iV^ (10) 

extends generally up to 1 -f 2 MHz. For the sake of 
simplicity it has been approximated by a uniform spec- 
trum extending up to 1 AlHz . 5v and 8T extend on very 
narrow bandwidths, while N Xo and w p are white noise 
sources. 

In Fig. r2Ja-d) we have plotted the five T L (9) vs. 9 
for three excitation strengths E = 0.71,0.87,0.975 of a 
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FIG. 2: Tg L (see Eq. (|D1[) ) vs. 9 for different noise sources: 
(a) pump quantum noise; (b) amplitude and temperature; (c) 
cavity length; (d) pump phase. The plots, referred to dif- 
ferent scaling factors, have been calculated for typical values 
of ft f (= 0.3), 7/ (= 0.15) and E = 0.71(dashed blue), 0.87 
(purple dotted), 0.975 (black). 



perfectly tuned OPO with a pump cavity mode linewidth 
twice the signal one (ko = 2), a condition similar to that 
of Ref. [2l| . The graphs show that for iV X0 , fi p , and ST the 
maximum deviation from a Gaussian appears for 9 = 0, 
while for w p and 8v it occurs for 9 = zkir/4. (see Eq. 




FIG. 3: Maxima of T 9t (see Eq. (|DT) vs. ko for the conditions 
of Fig. (E -> 0.975). 

The maxima of T L (9) for _/V Xo , w p and Sv decrease 
for increasing kq (see Fig. EJa-d)). As expected from 
Eq. (fT2"j) . the contribution of ST is independent of Kq. 
The same holds true approximately for jl p too having 
considered a technical noise bandwidth small compared 
to Ko- 

Looking at Fig. [3] we see that the maximum for t 
is at least one order of magnitude larger than the other 
ones. Moreover, g^ p ^> gx so that for pump level up to 
E 2 = 0.95, the NG behavior is essentially due to the laser 
amplitude fluctuations. 
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FIG. 4: T 0xo (pump quantum noise) and Tq Mp (amplitude 
fluctuations) vs. 1 — E 2 . 



In Fig. 1T 0X0 
are plotted vs. 1 - 



(blue-dashed) and T 0Mp (black-solid) 
E 2 on a double logarithmic scale. Ap- 



proaching the threshold the influence of iV Xo increases 
dramatically although To Xo does not overcome To^. 
Moreover, g x -C g^ p , so that the observation of pure 
quantum effects, predicted by Drummond et al. [26j, is 
demanded to future technology when either new mate- 
rials with huge non-linear coefficients (enhanced g x ) or 
very quiet lasers (reduced g^ p ) will be available. 

Being the laser amplitude noise the prominent source 
influencing the non-linear behavior we have compared 
some experimental findings of Ref. 2l| with the predic- 
tions of the herein discussed model. In particular, the 
experimental behaviors of Kg and Kg^ vs. 9 G (— tt, tt) 
are reported in Fig. [5] Moreover, the maximum value of 
the experimental kurtosis, for 9 = 0, and ff^ p To Mp (see 
Eq. (HH)) vs. E 2 e (0.45,0.97) are plotted in Fig. M 

In general T^^ can be represented by 



T^ p = T 4Mp cos 46* + T 2Mp cos 29 + T 0Mp 



(17) 



with T4 j 2 j o/j functions of E, ko and k. For very small 
deviations from the resonant configuration and \E\ close 
to 1, Y# Mp depends critically on \E\ 



\E\ 



(18) 



with Eq = e\x\/ [iplVj the excitation strength at res- 
onance. 

The different heights of the peaks in the experimental 
data (see Fig. [S]) can be can be ascribed to the variation 
of E during a measurement. The acquisition time for Vg 
for each 9 lasted about 2 ms, implying a total 9— scanning 
acquisition time of 200 ms. The apparatus was equipped 
with a digital controller providing a crystal temperature 
time constant > 10 3 s and a Drever-Pound system con- 
trolling the cavity tuning with a time constant > 10 s. 
However, during the total acquisition time, slow drifts of 
the average cavity detuning v p can occur thus inducing a 
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-tt/2 
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FIG. 5: K, 



(A) (see Ecl 
imental data of Ref. [2l|] 
72 - 0.92, n 



161) ) overimposed to the exper- 
Kfip (#) has been calculated by 
assuming E 2 = 0.92, Q f = 0.3, 7/ = 0.15 and g MP = 0.007. 
The horizontal error bar accounts for the detector phase 8 
stability, while the vertical one corresponds to the average 
spread between two neighbour Kg data. 



FIG. 6: Comparison of gf lp T 0llp (Eq. JT5J) vs. E 2 (solid 
curve; Qf = 0.3, 7/ = 0.15 and g Mp = 0.007) with exper- 
imental data. Each experimental point has been obtained 
multiplying the measured kurtosis at 6 — for the relative 
squared variance. It is evident the agreement between data 
and theoretical curve within error bars. 



variation of the effective excitation parameter E during 
the scan (see Eq. (|18p). The simplest and more direct 
way to account for the time dependence of v p is to set in 

v p = a{0- 0o) (19) 

with a and 9q two fitting parameters and express the 
coefficients T^.o^p m Eq. (|T7f as functions of 9. 

The best fit (continuous curve of Fig. [5]) of the ex- 
perimental data has been computed by assuming a spec- 
tral density uniform in the interval 4- 1.0 MHz, in 
agreement with the laser (Lightwave mod. 142) technical 
noise specification and optimizing Eq (see Q18p). a and 
6*o (see (HU)). The best agreement has been obtained for 
a = 0.013, 9q = 7r and E — 0.932, values corresponding 
to a drift of the resonance frequency of ~ 4% the pump 
cavity mode linewidth and a variation of \E\ of 0.006. Fi- 
nally, normalizing Tg^ to the squared experimental 
variance (see Eq. (fT6|l ) the best agreement was obtained 
for g^ = 0.007 in agreement with laser noise specifica- 
tion (< 1%). 

The NG character depends critically on the distance 
from the threshold (sec Fig. @J. gfc Tq Mp (see Eq. ([H)j)) 
vs. E 2 is compared to a set of 35 data obtained for five 
different values of E 2 (0.5, 0.7, 0.8, 0.9, and 0.95) in 
Fig. [5] Plotted values have been obtained multiplying 
the experimental kurtosis by the relative squared vari- 
ance. Error bars have been obtained by standard error 
propagation. The good agreement between the expected 
behavior and the data confirms, once more, the effective- 
ness of the model and the validity of assumptions about 
the relative noise strengths. 



V. CONCLUSIONS 

The statistical properties of the fields generated by an 
OPO depend on the fluctuations of many classical param- 
eters, namely pump amplitude (g Mp /tp), pump frequency 
(g-cvptVp), cavity detuning (g v 8v), and crystal tempera- 
ture {gxST). In this paper it has been presented a model 
of the OPO based on an extension of the Graham-Haken 
quantum Langevin equations (GHLE) which accounts for 
the fluctuations of these parameters. The field generated 
has been dealt with as the response of a nonlinear device 
to these noise sources. Then, expanding the extended 
GHLE system at different orders in g^ p , g VJp ,g u ,gT a hier- 
archy of equations has been obtained with p, p , vj p , Si), ST 
acting as noise sources together with the quantum noises 
(g x N x , g x N Xo ) entering the optical cavity. These sources 
have been modeled as Gaussian processes with unit 
standard deviations weighted by the respective factors 
gn P ,9vj P ,gv,9T typically ranging in the intervals 10~ 2 
10~\ IO- 7 / 2 ~ 10~ 2 , 10~ 5 -j- 10" 1 , 10~ 5 + 10~ 4 . 

The extended GHLE solutions, obtained beyond the 
linear approximation, have been used for assessing 
the non-Gaussian character of the field outing a de- 
generate OPO. The departure of the output quadra- 
ture X§ from the Gaussian statistics has been esti- 
mated by means of the Kurtosis- excess figure Kg = 

(( x t) ~ 3 (Xf ) 2 ) / (X|) 2 , i.e. the relative deviation of 
Xg 4-th moment from the Gaussian expression of it. The 
model furnishes Kg as a sum of contributions from differ- 
ent parameters. Sets of plots have been provided, show- 
ing the dependence of Kg on the OPO operating con- 
dition, namely, the ratio pump/signal bandwidths, the 
excitation strength E, and the detection frequency Q f 
and bandwidth 7/. For typical operating conditions the 
pump technical noise emerges as the most critical factor. 
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The model has been used for reproducing the exper- 
imental values, reported in Ref. [21] , of Kg vs. 9 and 
Kq vs. E 2 . The good agreement, within the error bars, 
of the experimental data with the analythic predictions 
confirms the validity of the presented model. By provid- 
ing a physically-insightful and computationally-effective 
parameterization of the OPO, the model may help in ad- 
dressing the generation of non-Gaussian states by means 
of OPO sources. 



Appendix A: Linearization below threshold 



In particular for the tuned case 



OTN (u) = 4 



l + E 2 



2E 
1 



2E 



(B2) 



Appendix C: B' 1 ' expansion 

The different noise sources N L contribute to B^ (Eq. 
T2|) ) through the terms: 

B« = g x N X0 B W + 9x NlM )T + E > ( C1 ) 



The Fourier transform of the Green's functions G, G 



relative to a degenerate OPO are given by 



G 



1 

D 



At 



\E\ 



Go 



V o 

o Ajr 1 



(Al) 

with $ = ip — ±Vo, A = k - iw, D (uj) = AA* - 
\E\ 2 = — (uj + u> + ) (uj + ujJ), A = kq - iw,uj± — 



i ^cos^T y\E\ —sin tp j . In particular for k = 1 
(tuned device) they correspond in the time domain to 




cosh (Et) sinh (Et) 
sinh (Et) cosh (Et) 



T > 







T < 



1 
1 



T > 

r < 



(A2) 



Appendix B: Time- normal ordered correlation 
matrix : cr' 1 ' 1 -* (ui) : 

The Fourier transform : c^ 1,1 ' (uj) : of the time-normal 
ordered matrix : a (r — r') :=: ^a' 1 ) (r) aW T ( T ')) x : 
(: a(r'-r) :=: a(r-r') : and : a Q t Q t :=: a* :, 
: a aa f :=: a a t a :) is obtained from cr^ 1,1 -' (a;) = 
a (w) / [(w 2 - uj\) (uj 2 - uj 2 _)] with 



a (uj) = 4 



A } A*(-l^)A(lj) 



1 



|B| 2 
A 



by first normally ordering a (uj) . 



a (uj) — > &n (uj) = 4 



A* 



1 

(A*)* 



and then, symmetrizing with respect to time reversal 



r(l,l) 



ctttv (w) 



2uj a 



1 2 



CT^v (-CJ+) - 



2w_ 



with 



A^ 1 




uj p 



e^\E\klk\- x 



e~ \E\ A- 1 



(i) 



B 



B« 



e~ M \E\ 
e^\E\ 



-1 



e~ M \E\ A n 



-e**|£?| Aj" 1 1 



(C2) 



Appendix D: Kurtosis— excess expansion 

From the vanishing of the time-normal ordered corre- 
lations (: a«Nf :) = (: a^Nj :) = 0, it follows that: 

(: X^X^ :) = e T - (: a^a^ T :) -6 

with I, m — 1,2. Hence, retaining only the lowest non 
linear orders for (: V e 4 :) and (: V g 2 :) 2 , Tg t (see Eq. ([To]) ) 
reads: 



-CTJV (- 



^ J-oo 

= T 4i .cos40 + T 2t cos26» + To t (Dl) 
with S L (w) = (^Ni (—w) Ni (w)\ and q L a 2 x 2 matrix 

1 1 f°° - 
= — / H (w, uj) a t (w, uj) du ,(D2) 

27T K - IW J_ 00 

T being a matrix transposition operator, <r L (w, lu) an 
entire function of u> and it) 

U-) 

cr L (w, lu) = (l + T^j (ko — iw) D (w + uj) G (w + uj) 
( B1 ) • BW (w) ■ 5tn (uj) 
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with cttn (w) defined in Eq. (IB1|) . 

H(wu)= Fj{w + w)F } {~u>) 

D{w + uj)(uj 2 - u%) (w 2 - u?_ ) 

and 

the Fourier transform of Ff (see Eq. (fl"5j) ) while fl± = 

Since lim cj i? (to , (J) cr L (w , to) — the RHS of Eq. 
(ID2[) is given by the sum of residues 

4 

&( w ) = _L_y;£rW(Te)* t («; >w ,) (D3) 
u i=l 

where (w) = Res u=u]l [H (w,uj)] for = 

cj + , w_, — n + ,— 0_ (J = 1, 2, 3, 4), poles of H (u>,w) in the 
upper complex w-plane. 

In the limiting case of zero centered delta-like sources 
Ni (see Eq. dUI) ) 

t Si = 9 T ■ q (o) • ee T ■ q (o) • e 

Such an approximation holds true for NttN v and in a 
less measure for N fip , depending on the laser technical 



noise bandwidth normalized to the OPO cavity one. On 
the other extreme, N Xo and N m represent white noises 
processes for which Yg XOtZa reduce to 



Tflxo,~,=-I> {,) H*°) 
i=i 

T • <r X a,™ P {-i&o,ui) ■ 0e T ■ S Xo ,a p (ik ) ■ 

4 1 r i 

- V m Res c0 H {1) (to) 

with the frequencies of (ID3|) and ty^ the poles of 
i?^^ (w) in the upper complex iz;-plane, 

= {2uj + , uj + — f2+, o; + — f2_, o; + +w_} 
= {-20+, 

Wl (4) = {-20_, lj_ - 0_, w+ - - fi_ - 0+} . 
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